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Production of open and hidden charm hadrons in heavy ion collisions is considered 
within the statistical coalescence model. The charmed quark-antiquark pairs are assumed 
to be created at the initial stage of the reaction in hard parton collisions. The number 
of these pairs is conserved during the evolution of the system. At the hadronization 
Q\ • stage, the charmed (anti)quarks are distributed among open and hidden charm hadrons 
in accordance with the laws of statistical mechanics. 

The model is in agreement with the experimental data on J/psi to Drell-Yan ratio in 
Pb+Pb collisions at SPS. This agreement can be reached only if a rather strong enhance- 
ment of the open charm production in central Pb+Pb collisions is assumed. A possible 
o ' mechanism of the charm enhancement is discussed. 

At the top RHIC energy, the model predicts an increase of J/psi to Drell-Yan ratio in 
more central nucleus-nucleus collisions with respect to less central ones. 
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The interest to charmed quarks in heavy ion physics was initially motivated by the 
suggestion of Matsui and Satz [p]] to use charmonia as a probe of the state of matter 
created at the early stage of the collision. 

The standard picture of charmonium production in nucleus-nucleus collisions assumes 
that charmonia are created exclusively at the initial stage of the reaction in primary 
nucleon-nucleon collisions. During the subsequent evolution of the system, the number of 
hidden charm mesons is reduced because of (a) absorption of pre-resonance charmonium 
states in the nuclei (normal nuclear suppression), (b) interactions of charmonia with sec- 
ondary hadrons (comovers), (c) dissociation of cc bound states in the deconfined medium. 
It was found that the J/ip suppression with respect to Drell-Yan muon pairs measured 
in proton-nucleus and nucleus-nucleus collisions with light projectiles can be explained 
by the normal nuclear suppression alone 0. In contrast, the NA50 experiment with a 
heavy projectile and target (lead- lead) revealed essentially stronger J/ip suppression for 
central collisions [^J. This anomalous J/ip suppression was attributed to formation of 
quark-gluon plasma [|J. 

Despite of quite successful agreement with the J/ip data, the standard scenario seems 
to be in trouble explaining the ip' yield. The recent lattice simulations || suggest that the 
temperature of ip' dissociation Td(ip') lies far below the deconfinement point T c : Td(ip') ~ 
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0.1-0.2T C 0. Therefore, not only the quark-gluon plasma, but also a hadronic co-mover 
medium should completely eliminate if)' charmonia in central Pb+Pb collisions at SPS. 
However, the experiment revealed a sizable if)' yield (see, for instance, [0]). It was observed 
U that if)' to J /if) ratio decreases with centrality only in peripheral lead-lead collisions, 
but remains approximately constant at sufficiently large number of participants N p > 100. 

A completely different picture of charmonium production was proposed in Ref. M: 
hidden charm mesons are supposed to be created at the hadronization stage. Similar 
to all other hadrons, their abundancies can be described within the thermal model [|T0fl . 
However, production of heavy quarks in soft processes is expected to be negligible. Most 
likely, they are produced exclusively at the hard stage. Their number, therefore, can, 
generally speaking, deviate from the thermal equilibrium value. This fact was taken into 



account by the statistical coalescence model (SCM), which was proposed in Ref.[|TT|]. Soon, 
it was realised that, when the average number of heavy quark pairs per collision is small 
(< 1), exact conservation of this pairs in the hadronization process appears to be crucial 



In this talk, I present the results [)13T , |14f obtained within the statistical coalescence 
model: 

• c and c are created at the initial stage of the reaction in primary hard parton 
collisions; 

• their number remains approximately unchanged during the subsequent evolution; 

• they are distributed over open charm hadrons and charmonia at the hadronization 
stage in accordance with laws of statistical mechanics. 

Within its applicability domain (N p > 100), the model demonstrates excellent agree- 
ment (x 2 /dof = 1.06) with the NA50 data on J/ifi production (see Fig. [[]). 

The so-called 'second threshold of anomalous J/if> suppression' at large transverse en- 
ergy (E T > 100 GeV) naturally appears in SCM due to two effects: (i) fluctuations of 
the transverse energy at fixed number of participant nucleons and (ii) -fir-losses in the 
dimuon event sample with respect to the minimum bias one (see Ref. |14| for details). 

The agreement with the data is reached assuming a rather strong enhancement of the 
open charm production in central Pb+Pb collisions: up to a factor of about 3.5 within 
the rapidity window of the NA50 spectrometer. This is consistent with the indirect 
experimental result [T5"|. 



A possible mechanism of open charm enhancement was proposed in Ref. fllEf . Calcu- 
lations in the leading order of perturbative quantum chromodynamics (QCD) show that, 
assuming factorization (i.e. independence of the hard QCD process on the subsequent 
soft stage), a great fraction of cc pairs is created with invariant masses M C c below the cor- 
responding meson threshold 2m £>. In vacuum, these subthreshold pairs would not be able 
to form observable hadrons. In reality, however, some part of these pairs may hadronize 
due to interaction with spectator partons. Nevertheless, there is no reason to assert that 
all subthreshold pairs result in final state hadrons with open or hidden charm. It is quite 
reasonable to assume that a great part of them has to annihilate into lighter hadrons, or, 
alternatively, the final state interactions prevent formation of such pairs (factorization is 
broken) . 
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Figure 1. The dependence of the J /if) to 
Drell-Yan ratio on the transverse energy at 
SPS. The vertical line shows the boundary 
of the applicability domain of the statisti- 



cal coalescence model (SCM), see [[LJ] for 
details. 
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Figure 2. The dependence of the J/ if) to 
open charm ratio on the number of nu- 
cleon participants at RHIC. The vertical 
line shows the boundary of the applica- 
bility domain of the statistical coalescence 
model (SCM). 



In presence of a deconfmed medium (quark-gluon plasma or its precursor), the sub- 
threshold quark pairs can hadronize due to interactions with the medium, or, alterna- 
tively, the medium can modify the final state interactions making possible creation of 
the subthreshold pairs and their subsequent hadronization. This should lead to enhanced 
production of the total charm in nucleus-nucleus collisions in comparison to the standard 
result obtained within the direct extrapolation of nucleon-nucleon data. 

We have estimated the upper bound of the open charm enhancement. At SPS energies, 
one can expect maximal enhancement by a factor of about 6, which is consistent with our 
above fit for the J /if) to Drell-Yan ratio. 

It is interesting to note, that the statistical coalescence model does not always lead to 
J /if) suppression. Predictions of the statistical coalescence model for RHIC are shown in 
Fig. [I] (see also Ref . [p]j ) . As is seen, J /if) suppression is still expected at the lowest RHIC 
energy. In contrast, relative J /if) enhancement (e.g. growth of the J /if) to open charm[] 
ratio with the number of participants at N p > 100) is predicted for the top RHIC energy. 

In conclusion, the NA50 data for not very peripheral lead-lead collisions are consistent 
with the following scenario: 



1 It should be stressed that, due to possible nuclear shadowing, the open charm multiplicity has not to be 
proportional to the number of binary nucleon-nucleon collisions. Therefore the ratio of J / if multiplicity 
to the number of binary collisions may decrease even if J/tf) to open charm ratio increases. 
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Formation of deconfmed medium leads to enhanced production of charmed hadrons. 
Charmonia as well as other hadrons are formed at the hadronization stage. The distribu- 
tion of charm quarks and antiquarks over open and hidden charm hadrons follows laws of 
statistical mechanics. 

A relative J/ip enhancement is expected at the top RHIC energy. 

I am indebted to M. Gorenstein, W. Greiner, L. McLerran and H. Stocker for fruit- 
ful collaboration. I wish to thank A. Andronic, F. Becattini, P. Bordalo, P. Braun- 
Munzinger, L. Bravina, A. Capella, Yu. Dokshitzer, J. Cleymans, W. Florkowski, A. Fraw- 
ley, L. Grandchamp, A. Kaidalov, B. Kampfer, J. Kapusta, F. Karsch, D. Kharzeev, 
V. Koch, K. Redlich, D. Rischke, M. Rosati, H. Satz, E. Shuryak, Yu. Sinyukov, J. Stachel, 
J. Wambach, A. Zhitnitsky and G. Zinovjev for interesting discussions. 
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